Gravitational dependence of contact angle hysteresis 
for a solid defect pinning a liquid-vapor interface 

J0rgen Vitting Andersen and Henrik Jeldtoft Jensen, 
Department of Mathematics, Imperial College, Huxley Building, 
180 Queen's Gate, London SW7 2BZ, United Kingdom 
(Last revised February 1, 2008) 



On 
OS 

o 

Q 



> 
o 

(N 
-I— > 



c 

o 
o 



X 



We study the case of contact angle hysteresis, when a solid defect is quasi-statically pulled out 
and pushed into a liquid, in the presence of gravity. We solve the problem analytically in an 
approximation which, in the case of a truncated parabolic pinning force, becomes exact in the the 
limit of slow variations in the height of the liquid vapour interface. The case of a Gaussian pinning 
force is solved by a geometrical construction. We show that the extension of the defect, and not the 
capillary length, is the determining length scale, when the extension of the defect is much smaller 
than the capillary length. 



PACS numbers: 



.lO.Cr, 68.45.Gd 



The spreading of a liquid on a solid is known to 
be strongly influenced by chemical inhomogeneities and 
roughness of the solid surface Different situations 
arise for the motion of the co-existence line between a 
solid, a liquid and a gas phase (called the triple line or 
the contact line) depending on whether the liquid com- 
pletely wets the solid surface m or incompletely wets 
the solid [3-18]. In the case of incomplete wetting, which 
we will consider in this letter, the onset of hysteresis for 
single defect without gravity was first solved analytically 
in ||,D. Much effort, experimentally [3-10] as well as 
theoretically [11-19]. has since then been done to un- 
derstand the case with gravity as well as the case where 
many defects pin the triple line. 

Here we study contact angle hysteresis for a single 
Gaussian defect in the presence of gravity. Expanding 
the Gaussian pinning force to second order, we derive 
an analytic expression for the contact angle hysteresis. 
In particular we find a weak logarithmic dependence on 
gravity when the extension of the defect is small com- 
pared to the capillary length. In this limit the hysteresis 
is mainly determined by the ratio between the liquid- 
vapor surface energy, and the strength of the pinning 
potential times the extension of the defect. Generalising 
to the case where the defect is placed on a solid plate, we 
show that the hysteresis then has an additional depen- 
dence of the capillary length because of the capillary rise 
along the solid plate. 

We use a macroscopic Hamiltonian for the system con- 
sisting of a (horizontal) liquid and a solid vertical plate 
with a fiber of defects. |ll| Let the liquid-vapor inter- 
face, h(x,y), be chosen along the s-y-direction, with the 
solid being pulled vertically out of the liquid in the z— 
direction, see Fig. la. Minimizing the Hamiltonian leads 
to the force balance: 

pgh(x, y) - j LG V 2 h{x, y)(l + (Vh(x, y)) 2 )^ = 



A' p AyS(x)S(y) + A° p 6(x) (1) 

7lg is the liquid-gas surface energy, A p the strength of 
the pinning potential of the fiber defect and .A? is the 
strength of the pinning potential of the plate KG], p is 
the mass density of the liquid, g the gravity constant and 
Ay the width of the fiber of defects. Solving Eq. (|l|) in 
the limit (Vh) 2 < 1 gives: 
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2l pg the capillary length and Kq the zeroth- 

order modified Bessel function. 

We consider first the case without plate-liquid inter- 
action, A® = 0. In this case the problem is rotational 
symmetric around the axis of the fiber. The force per 
unit length (for a given width, Ay) needed to lift the liq- 
uid the distance h(0, Ay) along the fiber is determined 
by Eq. (§) as, 
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Dispite the fact that Eq. m) is obtained for a fiber of 
width Ay, the equation also describes the case of a defect 
of width Ay and zero extension in the z-dircction. This 
is because the part of the fiber which is above/below the 
triple line does not enter in the force balance of the triple 
line described by h{x, y). The pinning force is in balance 
with the gravitational force determined by the weight of 
the liquid meniscus F g — pgV . The volume V of the 
liquid meniscus is given by an integration over Eq. (|^) 
with the assumed A p from Eq. (||). One finds 



1 



F ° = tSt M", Ay) ee M(0, Ay) (4) 

where we introduced the effective elasticity constant k. 
Asumming the pinning force is Gaussian, with it's center 
z d above 2 = 0, 



F p (z) = A p Ay exp {-[(z - z d )/R} 2 ), 



(5) 



one can like B determine the critical amplitude A c v for 
onset of hysteresis by equating the maximum slope of 
F p (z) to the slope of F g (z): 
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Experimentally, hysteresis is defined as the difference be- 
tween advancing and receding work done on the contact 
line. We define accordingly the hysteresis H by, 



H 



F p (z t (z d ))dz d + F p {z t {z d ))dz d (7) 
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where Zt(z d ) describes the dependence of the position of 
the triple line as a function of the position of the defect 
z d . Notice H can be different from since F p is not nec- 
essarily a single valued function of z d . This is due to 
discontinuities in the function z t (z d ) when A p > A c p . As 
depicted in Fig. lb, z t (z d ) undergoes a jump from posi- 
tion Z\ to zi for decreasing z d and a jump from Z3 to Z4 
for increasing z d . The positions z\, z 3 are determined by 
the condition dF p (z) / dz = k and the positions z 2l 24 are 
determined by the condition kz = F p (z). The integrals 
in Eq. (^) can be transformed into integrals over zt by use 
of the Jacobiant dzt/dz d — (k — dF p /dz t )/k |22| and we 
can express the hysteresis in terms of the jump positions 
z n , n = 1, • • • ,4. 
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In the case of a Gaussian force Eq. (^), it is easy numeri- 
cally to find the jump positions z n by solving the geomet- 
ric construction sketched in Fig. lb. The resulting hys- 
teresis calculated from Eq. (||) is shown in Fig. 2 versus 
pr, where k* — A p \f2/(R^fe) is the maximal elasticity 
constant giving hysteresis. The inset of Fig. 2 shows the 
jump positions z n . Notice that the hysteresis diverges as 
l/k when p- — > 0, because z 3 — > while at the same 
time \z%\, \z 2 \, \z 4 \ — > 00. 

A closed analytic form for the hysteresis can be ob- 
tained if we replace the Gaussian pinning force in Eq. (|^) 
by a truncated parabola 



F P (z) = { 



if \z - z d \ < R 
otherwise 



(9) 



In this case Eq. (g) leads to a hysteresis given by 

1 3R 2 C C 2 R 3 



H( 1 ,pg,Ay,R) = A p Ay[—- 

+ T2 {c ' 2 + 4/R2 - 4C/i?)3/2] (10) 

where C = k/(A p Ay). The effect of the gravitational 
field enters through k in C . Since the zero order Bessel 
function Kq{x) depends logarithmically upon x for small 
x we conclude that when Ay <C l c the hysteresis only 
depends logarithmically upon the acceleration of gravity, 
i.e., a very weak dependence. In this case the extension 
of the defect (Ay, R) becomes the relevant length scale. 
Notice that this is in strong contrast to the common belief 
that the capillary length is the determining length scale 
in contact angle hysteresis. Contrary, when the defect is 
spatially extended (e.g. in the case of a wall with A® 7^ 0), 
we show below that the hysteresis depends on powers of 
l c , i.e. a strong dependence of gravity. 

Like the Gaussian case we have plotted the hysteresis 
for the parabolic force in Fig. 2 versus pr. One notices 
that the two different pinning forces give rise to same 
the l/k divergence for a weak elastic constant pr — ► 0, 
since in this limit the hysteresis is determined by the 
jump solution 23, which coincide for the two different 
shapes of pinning force. On the other hand for pr — > 
1 all jump solutions z n become relevant in determining 
the hysteresis, and the specific form of the pinning force 
matters. 

In order to illustrate our findings for a typical set of ex- 
perimental values, we have plotted in Fig. 3 the hysteresis 
versus gravity pg for different values of defect interaction 
range R. We have used a jlg corresponding to water- 
air, and assumed a contact angle hysteresis of the defect 
of 20°. We have furthermore assumed a fixed given as- 
pect ratio of the defect given by Rj Ay = 0.1. The curves 
(e) — (h) illustrate the very weak dependence of hysteresis 
on gravity when one is not in a parameter regime close to 
onset of hysteresis {pg w 1000,g/(cm 2 s 2 ) corresponds to 
the density of water) . For comparison note that a change 
in R by one decade makes a change in H by two decades. 
On the other hand when l c w Ay a pronounced depen- 
dence of gravity happens as illustrated by (a) — (d) . For a 
fixed pg one notices from the curves (6) and (d) (or from 
Eq. (|10|)) that H goes through a maximum as a function 
of R. The same trend happens for the Gaussian pinning 
force as can be seen from the curves (a) and (c) . 

Let us finally turn our attention to the general case 
with Ap ^ 0. We can repeat the procedure above. The 
gravitational restoring force for system of a plate with a 
fiber is given by 



F, 
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Ly is the length of the system in the y direction. The last 
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term on the right hand side of Eq. ( |11| ) is the contribution 
from the solid plate alone, whereas the negative term 
enters because the gravitational force on the fiber is now 
changed due to the additional capillary rise of the plate. 
The hysteresis in the general case with a plate plus fiber 
can be found analytically in the case of a pinning force 
given by Eq. (g): 

H(^,pg,Ay,R) = 

+ ^[(C 2 + 4/i? 2 + 4d/i? 2 - AC/Rf' 2 - {^f' 2 ] (12) 

where d = Ap c Tr/[V2K (^^)A' p Ay]. The linear de- 
pendence of C\ on the capillary length l c leads to a lead- 
ing dependence on gravity in the form (pg)^ 1 / 2 . 

Conclusion: We have studied the case of contact angle 
hysteresis, when a solid defect with a Gaussian pinning 
force is quasi-statically pulled out and pushed into a liq- 
uid, in the presence of gravity. Expanding the Gaussian 
pinning force to second order and assuming slow varia- 
tions in the height of the liquid vapour interface, we have 
derived an analytical expression for the contact angle hys- 
teresis. In strong contrast to the common belief, we have 
shown that the extension of the defect, and not the cap- 
illary length, is the determining length scale, when the 
extension of the defect is much smaller than the capillary 
lcnth. 
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FIG. 2. Hysteresis H versus -p- determined for a Gaussian defect force using Eq. (^) (thick solid line) and for a truncated 
parabolic defect force using Eq. ( |lo| ) (thin solid line). Inset shows jump positions z„ versus pr 
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FIG. 2. Inset to Fig. 2 
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FIG. 3. Hysteresis H versus pg for different values of interaction range R determined for a Gaussian defect force using Eq. (^) 
(thick solid lines) and for a truncated parabolic defect force using Eq. ( p^| ) (thin solid lines). We have used ~/lg = 73dyn/crn 
corresponding to water, and a typical experimental value, 20°, for contact angle hysteresis of the defect, giving A p = 69dyn/cm. 
We have used a fixed aspect ratio of the defect R/Ay = 0.1 (a-b) R = 0.02 cm, (c-d) R = 0.01 cm, (e-f) R — 0.001 cm, (g-h) 
R = 0.0001 cm. 
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